Cosmology with bulk pressure 
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Abstract 

It is well known that cosmological particle production processes may 
phcnomenologically be described in terms of effective viscous pressures. 
We investigate the consistency of this approach on the level of relativis- 
tic kinetic theory, using a Boltzmann equation with an additional source 
term that leads to a change in the number of gas particles. For a sim- 
ple creation rate model we find the conditions for collisional equilibrium 
of a Maxwell-Boltzmann gas with increasing particle number and discuss 
their cosmological implications. We also comment on the possibility of 
a bulk pressure driven inflationary phase within causal, nonequilibrium 
thermodynamics . 

1 Introduction 

Entropy producing processes have played an important role in the evolu- 
tion of the early universe. On the level of fluid cosmology, the simplest 
phenomenon connected with a nonvanishing entropy production is a bulk 
viscosity. A bulk viscous pressure is the only dissipative effect that is com- 
patible with the symmetry requirements of the homogeneous and isotropic 
Friedmann-Lemaitre-Robertson- Walker models. For the possible physical 
origin of bulk pressures in the expanding universe see |lj and references 
therein. 

On the other hand, there exist attempts to use effective viscous pressures 
in order to model particle production proceses in the early universe (see 
references in j^]). For certain quantum processes in the early universe 
'viscosity functions' have been calculated, opening the possibility of an 
effective imperfect fluid description of these phenomena. 
Generally, any production process will give rise to source terms in the 
macroscopic fluid balances. The interpretation of these contributions as 
effective bulk pressures uses the fact that a source term in the energy bal- 
ance of a fluid may formally be rewritten in terms of an effective bulk 
pressure. The advantage of a rewriting like this is obvious: An energy- 
momentum balance with a nonzero source term violates the integrability 



conditions of Einstein's equations, while the energy- momentum tensor of 
an imperfect fluid is a reasonable right-hand side of the field equations. 
Moreover, this approach allows us to calculate the backreaction of the cor- 
responding process on the cosmological dynamics. 

We shall investigate here whether the approach of regarding particle pro- 
duction processes as equivalent to viscous pressures is consistent with a 
kinetic description. In section 2 we sketch the general kinetic theory of a 
simple gas with nonconserved particle number and introduce an effective 
rate approximation for the production term. The conditions for collisional 
equilibrium are found in section 3, while the above mentioned effective 
viscous pressure approach is established in section 4. The backreaction 
of particle production processes on the cosmological dynamics under the 
conditions of collisional equilibrium is considered in section 5. In section 
6 we drop the equilibrium conditions and comment on the possibility of 
a bulk pressure driven inflationary phase within causal, nonequilibrium 
thermodynamics . 



2 Basic kinetic theory 

The one-particle distribution function / = / (x,p) of a relativistic gas with 
varying particle number is supposed to obey the equation Q 

l [/] ee pV,< -itpV |£ = c [f] + s (x, P ) . (i) 

L [/] is the Liouville operator and C[f] is Boltzmann's collision term. The 
source term S takes into account the fact that the distribution function / 
may additionally vary due to creation or decay processes, supposedly of 
quantum origin. On the level of classical kinetic theory we regard this term 
as a given input quantity. Assuming a linear coupling of the source term 
S to / and requiring a macroscopic description of the production process 
in terms of the zeroth, first and second moments of / only, provides us 
with the expression 

S(x,p)=(-^+v(x) S jf(x ) p) ) (2) 

for the source term in (|l|), where u a is the macroscopic 4- velocity. The 
creation process is characterized by the spacetime functions r and v. 



3 Collisional equilibrium 

As in the case 5* = 0, the condition for collisional equilibrium implies the 
structure 

/ (Z, p) = f (x, p) = exp [a + p aP a ] , (3) 

of the distribution function. Using the latter expression together with (g) 
in (jl|) yields the equilibrium conditions 

a = - , hfa a = 0, (3 {a , b) = — ^g a b , (4) 

which are less restrictive than in the case of particle number conservation. 
Different from the case without particle production (r^ 1 = v = 0) we find 
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that a may change along the fluid flow lines and that (3 a = u a /T obeys 
a conformal Killing equation for m ^ 0. The latter property implies the 
temperature behaviour T/T = — 0/3 that coincides with the temperature 
law for massless particles (radiation) (9 = u a a is the fluid expansion). 
We conclude that radiation and matter in the expanding Universe may be 
in equilibrium at the same temperature, provided the particle number of 
the matter component is allowed to change at a specific rate. It is well 
known, that for conserved particle numbers an equilibrium between both 
components is impossible. 

A further consequence of the equilibrium conditions (Q) is a modification 
of Tolman's relation [^). For nonvanishing particle production the quantity 
T^—goo is allowed to be time dependent. However, like a in (Q), it has to 
be spatially constant. 



4 The effective viscous pressure approach 

Let, as usual, the energy-momentum tensor T %k be defined as the second 
moment of the distribution function (^). The effective viscous pressure 
approach then amounts to the replacement of the nonconserved energy- 
momentum tensor T lk = pu l u k + ph lk with p + (p + p) = p/r + vn by 
the conserved energy-momentum tensor T %k — pu % u k + (p + n)h lk with 
p + (p + p + 7r) =0, where 07r = — p/r — vn. We find this approach to 
be consistent for homogeneous spacetimes but not necessarily for inhomo- 
geneous ones M. 



5 Backreaction on the cosmological dynam- 
ics 

A nonvanishing effective bulk pressure influences the entire cosmological 
dynamics. Restricting ourselves to the homogeneous, isotropic and spa- 
tially flat case, the relation 2H / H = p/p is valid, 

where H is the Hubble parameter H = 0/3. The functions v and 
r are fixed by the conditions (4) and by the additional assumption of 
'adiabatic' particle production [3]. With H = a/a, where a is the scale 
factor of the Robertson- Walker metric, we find (see ||) that for massive 
particles (m ^ T) the scale factor behaves like a cx t 4 ^ 3 instead of the 
familiar a oc t 2 / 3, for r _1 = v = 0. Massive particles with adiabatically 
changing particle number are allowed to be in collisional equilibrium only 
in a power-law inflationary universe. Introducing a '7-law' in the usual 
form p = (7 — l)p, one obtains accelerated expansion, i.e., a/a > 0, in the 
interval 1 < 7 < 7 cr with 7„. w 1.09. In the opposite limit of massless 
particles the equilibrium conditions require r _1 = v = 0, i.e., there is no 
adiabatic production of relativistic particles. 

We also investigated whether a phase of exponential inflation is consistent 
with the conditions for collisional equilibrium In our setting expo- 
nential inflation may be realized by 'nonadiabatic' particle production. It 
turns out that the equilibrium conditions ^) are compatible with a de Sitter 
phase only for a time interval of the order of the Hubble time. For larger 
times exponential inflation violates the second law of thermodynamics, 
i.e., the de Sitter phase becomes thermodynamically unstable. 
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6 Causal thermodynamics and inflation 



Up to now our analysis was based on the equilibrium conditions (Q). Now 
we drop these conditions and investigate the role of bulk pressures within 
causal, nonequilibrium thermodynamics (see |^|, [sj and references therein) . 
The viscous pressure becomes a dynamical degree of freedom on its own 
and the cosmological dynamics is determined by a second-order equation 
for the Hubble parameter H . A basic ingredient of the theory is a nonva- 
nishing relaxation time. There exists a debate on whether deviations from 
thermodynamical equilibrium may drive a de Sitter phase. Assuming the 
applicability of the causal, nonequilibrium thermodynamics also far from 
equilibrium, it may be shown that the causal evolution equation for H 
admits stable, inflationary solutions with relaxation times of the order of 
the Hubble time, i.e., the nonequilibrium may be 'frozen in'. Provided, 
the bulk pressure is again due to particle production processes, one finds 
that the temperature remains approximately constant during the de Sitter 
phase (no supercooling!). The latter feature is connected with an exponen- 
tial increase of the comoving entropy. While the standard scenario of an 
inflationary universe comprises an adiabatic de Sitter phase and a subse- 
quent reheating period, in which all the entropy in the universe is produced, 
the entropy production takes place during the de Sitter stage itself in the 
present picture. 
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